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The  enclosed  reprints  and  preprints  represent  work  carried  out  by  personnel  sup¬ 
ported  under  this  contract.  Summaries  of  this  work  along  with  other  work  performed 
but  not  as  yet,  published  are  contained  within  this  report. 

We  (Stegeman  and  Seaton)  were  invited  to  produce  a  review  paper  on  Degenerate 
Four-Wave  Mixing1  with  guided  waves  which  was  published  in  a  special  issue  of  IEEE 
Journal  of  Quantum  Electronics.  One  of  the  conclusions  of  that  review  was  the  need 
development  and  understanding  of  new  nonlinear  optical  materials  suitable  for  nonlinear 
guided  wave  devices. 

In  light  of  this  need  a  major  part  of  our  research  program  in  the  last  year  has  been 
the  continuation  of  research  into  the  semiconductor-doped  glass  material  system.  Bulk 
DFWM  and  photoluminescence2  experiments  on  these  materials  determined  that  they 
have  a  large  nonlinear  index  of  refraction  (n2  ~  10“u  m2/w)  with  fast  response  times 
(r  <  16  psec)  required  for  high  speed  optical  signal  processing.  Along  with  collaborators 
in  the  University  of  Glasgow  we  successfully  developed  a  successful  ion-exchange 
(Na+:K+)  technique  for  producing  low-loss  optical  waveguides  in  these  materials,  the 
details  of  which  are  in  the  enclosed  reprint.2 

We  have  been  studying  the  effects  of  incorporating  "real"  material  parameters  into 
the  operation  of  nonlinear  optical  waveguide  devices4  and  have  determined  that  the 
most  important  parameter  for  guided  wave  devices  is  AnMt  ie.  the  maximum  change  in 
refractive  index  which  can  be  induced  in  a  material  before  saturation  effects  set  in.  To 
determine  this  parameter  (which  has  not  previously  been  studied)  we  have  been  per¬ 
forming  bulk  DFWM  experiments  in  a  number  of  semiconductor-doped  glass  samples. 
By  plotting  DFWM  reflectivity  against  pump  intensity  we  observe  signal  saturation  over 
a  range  of  samples  at  ~  1  MW/cm2  intensity  levels.  Using  both  a  two-level  atom  satu¬ 
ration  model  and  on  exponential  decay  model  we  determined  Ansat  of  1.7  x  10'4  and 
1.13  x  10-4  respectively.  Experiments  are  underway  to  determine  the  origin  of  "per¬ 
manent"  gratings  which  appear  to  be  created  via  a  photo  darkening  effect  in  the  semi¬ 
conductor  glasses.  This  effect  appears  to  be  particularly  prevalent  in  the  Schott  soda- 
lime  glasses  which  we  specially  obtained  to  facilitate  the  ion-exchange  process.  The 
problem  of  "solarization"  ie.  photo-darkening  is  known  to  be  particularly  prevalent  in 
soda-lime  based  glasses  so  color  filter  samples  with  borosilicate  host  glasses  (ie.  BK.7  or 
standard  Corning  filters)  will  be  further  investigated  for  ion-exchangability. 


.  Other  work  which  has  been  performed  with  AFOSR  support  from  this  contract  has 
been  of  a  theoretical  nature  regarding  nonlinear  guided  wave  (NLGW)  structures  where 
the  nonlinearity  is  sufficiently  large  to  effect  the  guided  wave  field  distribution.  We 
have  demonstrated  numerically  the  efficient  excitation  of  stable  TEq  NLGW's  with  a 
suitably  chosen  gaussian  input  beam.  Selective  excitation  of  any  of  three  possible  field 
distributions  corresponding  to  the  same  flux  level  in  a  system  comprising  a  linear  thin 
film  waveguide  bounded  by  two  nonlinear,  self-focussing  media  was  demonstrated 
numerically.5 

We  also  observed  that  external  excitation  of  a  NLGW  can  produce  sequential  thres¬ 
hold  behaviour  through  multi-soliton  emission  from  the  nonlinear  wveguide.6  This 
effect  has  intriguing  consequences  for  original  optical  switching  devices. 
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Abstract— The  theory,  experimental  implementation,  and  applica¬ 
tions  of  degenerate  four-wave  mixing  with  guided  waves  are  reviewed. 

I.  Introduction 

DEGENERATE  four-wave  mixing  (DFWM)  has  been 
studied  extensively  in  bulk  media  for  a  number  of 
years,  and  many  applications  have  been  predicted  and 
demonstrated  [1].  In  the  most  general  interaction  geom¬ 
etry.  there  are  three  input  beams,  two  of  which  are  coun- 
terpropagating  and  a  third  signal  beam  incident  at  an  ar¬ 
bitrary  angle,  all  of  which  interact  to  generate  a  fourth 
wave  via  the  third-order  susceptibility  of  the  medium.  The 
process  is  automatically  phase  matched  and  the  conjugate 
beam  is  emitted  back  along  the  direction  of  the  incident 
signal  beam.  The  efficiency  of  the  process  is  proportional 
to  the  product  of  the  two  pump  beams,  the  length  of  the 
interaction  region  along  the  direction  of  the  conjugate 
beam,  and  the  magnitude  of  the  nonlinear  susceptibility. 

The  electromagnetic  waves  guided  by  the  interface  be¬ 
tween  two  semi-infinite  media,  by  single  or  multiple  films 
bounded  by  two  semi-infinite  media,  or  by  fibers  can  be 
used  to  implement  DFWM.  The  key  feature  introduced 
by  guided  waves  is  that  the  fields  decay  exponentially 
away  from  the  boundaries  into  the  semi-infinite  bounding 
media  with  1/e  distances  of  typically  a  fraction  of  the 
wavelength  of  the  radiation  being  guided.  Therefore,  the 
effective  beam  dimension  along  the  direction  normal  to 
the  surface  can  be  on  the  order  of  the  wavelength  of  light, 
which  corresponds  to  the  minimum  beam  cross-sectional 
area  and  hence  maximum  power  density  for  a  given  input 
power  level.  Furthermore,  these  high  fields  can  be  main¬ 
tained  for  long  propagation  distances,  centimeters  in  the 
case  of  thin-film  or  channel  waveguides,  and  meters  or 
kilometers  for  optical  fibers.  Thus,  the  high  fields  and  long 
phase-match  distances  required  for  efficient  interactions 
can  be  achieved  with  modest  or  low  total  powers.  These 
attractive  features  were  recognized  at  an  early  stage  and 
have,  by  now.  been  applied  to  an  impressive  number  of 
nonlinear  guided  wave  phenomena  [2]. 
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Guided  waves,  in  general,  result  from  a  coupling  be¬ 
tween  an  electromagnetic  wave  and  some  resonance.  The 
resonances  can  be  geometric,  for  example,  as  occurs  in  a 
fiber  (or  thin-film)  waveguide  [3],  [4}  where  constructive 
interferences  across  the  core  or  film  result  in  waveguide 
modes.  The  resonances  can  also  be  material  properties  re¬ 
lated.  For  example,  the  plasma  resonance  associated  with 
the  electron  gas  in  a  metal  is  coupled  to  an  electromag¬ 
netic  field  via  the  interaction  between  the  field  and  the 
charges  and  leads  to  surface  plasmons  [5].  Such  material 
resonances  are  usually  accompanied  by  losses  and  there 
is  always  a  tradeoff  between  the  propagation  distance  and 
the  field  confinement  [6).  For  this  reason,  most  work  in 
guided  wave  nonlinear  optics  for  device  applications  has 
been  limited  to  dielectric  waveguides  with  freely  propa¬ 
gating  guided  waves. 

The  first  experiments  reported  [7]  on  waveguide 
DFWM  were  performed  in  CS2  filled  fibers  in  which  the 
long  interaction  distance  was  more  than  enough  to  com¬ 
pensate  for  the  small  nonlinearity  associated  with  CS2. 
The  phase  conjugate  properties  of  the  interaction  were 
demonstrated.  Despite  promising  applications  [8]  such  as 
wide-angle  narrow-band  filtering,  image-frequency  con¬ 
version,  and  small-sample  spectroscopy,  there  has  been 
little  progress  in  phase  conjugation  in  fibers  since  about 
1980.  Recent  interest  has  centered  on  experiments  in 
planar  geometries;  specifically,  experiments  using  thin- 
film  guided  waves  [9]  and  surface  plasmons  [10]  were  re¬ 
ported.  Because  of  the  potential  applications  to  "real- 
time”  signal  processing  [2],  [10].  these  initial  experi¬ 
ments  represent  the  first  steps  in  a  field  which  should  de¬ 
velop  rapidly  over  the  coming  years. 

In  this  paper,  we  review  progress  in  degenerate  four- 
wave  mixing  with  guided  waves.  Although  the  emphasis 
is  on  experimental  results,  the  basic  formulation  of  the 
problem  and  its  applications  will  also  be  discussed. 

II.  Theoretical  Formulation 

Degenerate  four-wave  mixing  can  be  carried  out  with 
either  two  or  four  beams  guided.  Certainly,  the  most  ef¬ 
ficient  interactions  occur  with  all  four  waves  guided  since 
this  configuration  optimizes  the  power  density  in  each 
wave  for  a  given  power.  For  the  two-wave  guided  case, 
two  externally  incident  waves  can  be  used  to  phase  con¬ 
jugate  a  guided  wave  [12].  or  two  guided  waves  can  be 
used  to  phase  conjugate  an  externally  incident  plane  wave  I 
field  [13].  Here  we  consider  only  the  four  guided  wave 
case. 
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A.  Guided  Waves 

Guided  waves  are  electromagnetic  fields  which  satisfy 
both  Maxwell’s  equations  in  every  medium  into  which  the 
fields  penetrate  and  continuity  of  the  tangential  boundary 
conditions  at  every  interface.  For  the  simplest  case  of  iso¬ 
tropic  waveguide  media  with  propagation  along  the  x  axis, 
the  sth  input  wave  (s  =  1,  2,  3)  in  the  mth  normal  mode 
an  be  written  as 

£(r,  t)  =  4  E{m  s>(r)  a'n  s'(x) 

■  exp  [/(a>r  -  pm  s%  ■  r)]  +  c.c.  (1) 

where  0im-s)  is  the  effective  index  obtained  from  the  so¬ 
lution  of  the  appropriate  dispersion  relation  and  Ems)(r) 
is  the  corresponding  normal  mode  field.  Since  these  so¬ 
lutions  are  normal  modes,  they  satisfy  the  orthogonality 
relations 

f  dr+  E{m\r)  •  £w‘*(r)  =  „  &m  n  (2) 

J  —do 

where  the  vector  rt  is  orthogonal  to  the  propagation  direc¬ 
tion  defined  by  k  and  Am  n  is  the  normalization  constant 
chosen  so  that  |a("I)(x)|2  is  the  guided  wave  power  for  two- 
dimensional  waveguides  or  the  guided  wave  power  per 
meter  of  wavefront  for  planar  guided  waves.  These 
expressions  are  valid  for  all  guided  wave  fields,  fibers, 
channel  waveguides,  surface  plasmon  waves,  etc. 

The  mixing  of  three  input  guided  waves  leads  to  non- 
near  polarization  fields  which  are  usually  expanded  in 
products  of  the  mixing  fields.  A  nonlinear  polarization 
field  of  frequency  o»  is  generated  in  each  of  the  media  in 
which  the  guided  wave  fields  exist.  It  is  usually  written 
as 


P(NrL„  =  |  P$  exp  [i(wr  -  0pkpx)]  +  c.c.  (3) 


where  (3pkp  is  the  polarization  wavevector.  The  generation 
of  the  nonlinearly  produced  conjugate  wave  can  be  treated 
’  •  the  slowly  varying  amplitude  approximation  using  well- 
lown  guided  wave  coupled  mode  theory  [3],  [4]  with  the 
result  that,  for  propagation  along  the  x  axis. 


(x)  =  i 


l$n-s)Ain 


dr,  ■  E*{n  s\r) 


•  exp  [~i(0p  -  0{n  s})x\.  (4) 

For  a  fiber  or  a  channel  waveguide,  the  integral  is  taken 
over  the  cross-sectional  area,  and  in  the  case  of  planar 
waveguides,  the  integral  is  performed  ove.  the  depth  co- 
"dinate. 


Fig.  1.  The  planar  four-wave  mixing  geometry  discussed  in  this  paper. 
The  pump  beams  are  I  and  2.  the  signal  wave  is  3.  and  the  conjugate 
wave  is  4  [14]. 


with  wavevector  as  shown  in  Fig.  1  for  a  planar 

guided  wave  geometry.  Wavevector  conservation  is  given 
by 

+  0(m'-2)k2  +  0,m'-3>*3  +  =  0  (5) 

with  k  |  =  k2  =  kj  =  kA.  If  the  two  input  beams  1  and  2 
are  exactly  contradirectional  and  have  the  same  mode 
number  m  =  m\  0im  ])k{  =  -&m-2,k2.  Therefore,  it  is 
necessary  that  m"  =  n  so  that  0'’lA)k4  =  -0ln-3,ki,  and 
signal  beam  4  is  generated  backwards  along  the  direction 
of  incidence  of  beam  3.  independent  of  the  incidence  di¬ 
rection  of  beam  3.  (For  a  fiber,  all  four  beams  propagate 
along  the  fiber  axis.)  Note  that  this  process  does  require 
that  beams  1  and  2  have  the  same  mode  order  m,  that 
beam  4  is  radiated  into  the  same  mode  index  n  as  beam 
3,  and  m  and  n  are  not  necessarily  equal  [14],  Because 
there  are  multiple  guided  wave  solutions  associated  with 
the  same  guided  wave  polarization  characteristics,  the 
spectrum  of  possible  interactions  is  larger  than  that  for 
plane  waves.  (The  other  consequence  is  that  the  possible 
directions  of  the  incident  waves  are  restricted  to  one  (fiber 
or  channel  waveguide)  or  two  (planar  waveguide)  dimen¬ 
sions  for  an  all-guided  wave  interaction.) 

For  an  isotropic  material,  the  conjugate  polarization 
term  [14]  is  given  by 

P  (r)  =  [4f()\  nj’f- w.  u>.  —w.  u>) 

•  £*'"V>  E)m  2\r)  E\m  U{r) 

+  2«ox  r::i(  — «.  iji)  £jm  "(r) 

•  E)m2\r)  £*'"  V)1  a{(x)  a2(x)  af(x )  (6) 

where  the  notation  for  the  field  amplitudes  has  been  sim¬ 
plified  so  that  a^x)  ■  ain  "(x).  This  expression  simplifies 
far  from  any  resonant  behavior  in  the  xl3)  terms  and 
Xna  =  X  i22i  =  Xeff-  For  guided  wave  fields. 

PNL”(r)  =  2«5n;n:>  [$  V)  K  "(r)  E\m/2\r) 

+  \  E*(;  h{r)  £«--2'(r)  ££  "(r)] 


B  Degenerate  Four-Wave  Mixing  in  Guided  Wave 
Geometries 

The  same  formalism  can  be  used  to  discuss  the  four- 
wave  mixing  in  one-  and  two-dimensional  waveguides  [8], 
[14],  There  are  three  input  beams  with  wavevector 
{&"  [0*"  2)*2],  and  [&m  -3)ifcj]  and  one  output  beam 


•  a,(x)  a:(x)  af(x)  (7) 

in  each  medium  (labeled  7)  where  *  0.  Here  we  have 
defined  3xefr  =  where  n2 T  is  the  intensity-de¬ 

pendent  refractive  index  for  the  yth  medium.  Clearly, 
similar  equations  can  be  derived  for  the  nonlinear  polar¬ 
izations  driving  the  pump  and  signal  beams. 
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Here  we  consider  only  the  two  cases  which  have  been 
studied  experimentally.  For  fibers  or  channel  waveguides, 
the  coupled  mode  equations  are 

—  a4(x)  =  (<K| ,22  +  *12:1)  a*(*)  a2(x)  a,U)  (8a) 
ax 

d 

—  a}( x)  *  -i(*n22  +  *1221)  <*4(Jt)  aiU)  a2(x)  (8b) 

dx 
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In  the  weakly  guiding  case,  that  is,  one  in  which  the  index 
difference  across  the  core-cladding  interface  for  fibers  or 
across  the  channel -bounding  media  interfaces  is  small. 
*1221  -  *1122  And  (8)  simplify  to  the  well-known  equations 
for  four-wave  mixing  with  plane  waves  with,  of  course, 
*1122  given  by  the  “overlap  integrals.”  In  the  absence  of 
pump  depletion  with  the  boundary  conditions  u3(0)  =  a, 
and  a4(L)  =  0  where  L  is  the  length  of  the  interaction 
region,  the  results  are 


a2(x)  =  a,  ,  tan  —  sin  —  +  cos  — ; 
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Consider  also  the  case  of  waves  guided  by  planar  sur¬ 
faces.  Assuming  that  TE  beams  1  and  2  propagate  along 
x  ’  which  is  oriented  at  an  angle  9  to  the  x  axis  along  which 
TE  beams  3  and  4  are  traveling, 

£W.ior2)(r,)  m  E,n  '\r,)  [sin  6.  cos  8.  0J; 

£i*.Jor«»(r)  =  £<*.3 nr)  f0  !  oj  (12) 


Again,  in  the  limit  of  no  attenuation  or  depletion  for  the 
input  beams,  and  assuming  that  beams  1.  2.  and  3  are 
introduced  into  the  interaction  region  at  x  ‘  =  0,  x'  =  L‘. 
and  x  =  L,  respectively,  substituting  into  (8),  and  inte¬ 
grating  the  signal  beam  from  0  —  L  gives  the  same  form 
as  (10)  and  (11),  but  with 
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In  the  general  case,  it  is  necessary  to  take  incident  beam 
depletion  into  account  which  leads  to  a  series  of  coupled 
mode  equations  between  the  amplitudes  of  the  various 
beams. 


III.  Experiments  in  Fibers 

In  1975,  Yariv  [15]  suggested  using  nonlinear  mixing 
as  a  means  of  compensating  for  phase  abberations  in¬ 
curred  in  transmitting  an  image  down  an  optical  fiber.  He 
subsequently  pointed  out  the  possibility  of  performing  this 
phase  conjugation  by  degenerate  four-wave  mixing  in  a 
multimode  waveguide  [16]  and  suggested  using  a  hollow 
core  fiber  filled  with  CS;  as  the  nonlinear  mixing  medium. 
Experiments  based  on  this  idea  were  reviewed  by  Hell- 
warth  in  1982  [17],  and  to  our  knowledge,  there  have  been 
no  further  results  published  since  then.  Thus,  the  only  ex¬ 
periments  to  have  been  done  are  those  of  Jensen  and  Hell- 
warth  [7]  and  Au  Yueng  ei  al.  [18]  in  CS2  filled  fibers. 
The  first  authors  achieved  high-fidelity  conjugation  with 
multimode  pump  waves  using  a  Q-switched  ruby  laser  and 
achieved  DFWM  reflection  efficiencies  of  =*50  percent 
for  35  kW  peak  pump  beam  powers.  They  also  demon¬ 
strated  that  the  wave  being  conjugated  could  serve  simul¬ 
taneously  as  its  own  pump  wave.  Au  Yueng  et  al.  [18] 
used  single-mode  pump  waves  to  obtain  0.45  percent 
DFWM  conversion  efficiency  for  only  6  mW  of  CW  pump 
power  from  an  argon  ion  laser.  However,  the  phase  con¬ 
jugate  nature  of  the  DFWM  signal  was  not  ascertained. 

In  related  experiments,  image  transmission  through  a 
multimode  optical  fiber  has  been  demonstrated  [19]  by 
using  a  photorefractive  BaTiO.  crystal  to  correct  for  modal 
distortions  external  to  the  fiber.  Stimulated  Raman  [19] 
and  Brillouin  [20]  scattering  in  fibers  can  also  be  used  to 
achieve  phase  conjugation,  but  these  topics  are  outside 
the  scope  of  this  paper. 

IV.  Experiments  with  Planar  Glided  Waves 

Two  experiments  on  four-wave  mixing  in  integrated 
optics  waveguides  have  been  reported  to  data  [9],  [21]. 
The  first  deals  with  optically  induced  damage  in  lithium 
niobate  waveguides  [21].  A  TE  guided  wave  was  created 
due  to  stray  scattering  via  the  photorefractive  effect,  and 
the  mixing  of  the  incident  TM  wave  with  this  TE  light 
creates  a  grating  which  further  enhances  the  TM  —  TE 
conversion.  Since  the  gratings  formed  are  written  more  or 
less  permanently  into  the  waveguide,  this  does  not  cor¬ 
respond  to  the  usual  four-wave  mixing  phenomenon.  The 
characteristic  time  scale  was  typically  seconds. 

The  second  experiment  was  performed  using  freely 
propagating  guided  waves  [9].  The  nonlinear  medium 
CS;.  which  was  used  as  a  cladding  on  top  of  a  thin-film 
waveguide,  was  located  inside  a  cell  optically  contacted 
to  the  surface  of  a  thin-film  glass  waveguide,  as  shown  in 
Fig.  2.  Three  coupling  pnsms  were  required  for  the  four 
beams  needed  in  the  interaction,  leading  to  very  difficult 
alignment  problems.  Since  the  nonlinear  mixing  occurred 
via  the  evanescent  tails  of  the  three  incident  guided  waves 
and  because  the  nonlinearity  of  CS:  is  small  (n;,.  »  3  x 
10' '*  m2/W).  the  reflectivity  (fractional  conversion  of 
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Fig.  2.  Sample  and  coupling  geometry  for  degenerate  four-wave  mixing 
in  a  thin-film  waveguide  showing  the  (hin-film  waveguide  (Coming  7059 
sputtered  onto  a  glass  microscope  slide),  the  arrangement  of  high-index 
coupling  prisms,  and  the  glass  cell  for  holding  the  CS,.  The  probe  beam 
intersects  ihe  pump  beams  at  90°  under  the  liquid  cell  (9). 


Fig.  3.  The  application  of  four-wave  mixing  to  (a)  convolution  (of  wave¬ 
forms  1  and  2)  and  (b)  time  inversion  lof  pulse  1)  [II]. 
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beam  3  into  beam  4)  was  only  =*  10-9.  Nevertheless,  the 
cubic  dependence  of  the  four-wave  mixing  signal  was  ver¬ 
ified  experimentally  for  15  ns  long  pulses  from  a  Q- 
s  witched,  frequency-doubled  Nd :  YAG  laser.  The  ob¬ 
served  signal  was  consistent  with  calculations  based  on 
the  reorientational  nonlinear  mechanism. 

V.  Planar  Waveguide  Applications 

This  interaction  has  some  possible  applications  to  real¬ 
time  all-optical  signal  processing  [11],  two  of  which,  con¬ 
volution  and  time  inversion,  are  illustrated  in  Fig.  3.  The 
two  input  signals  to  be  convolved  are  beams  1  and  2  which 
have  pulse  envelopes  (/,(r  -  x/v)  and  (A(f  +  x/v). 
Whether  the  convolution  process  occurs  or  not  can  be 
controlled  by  beam  3  which  is  assumed  to  have  a  constant 
amplitude  during  the  overlap  of  pulses  1  and  2  inside  the 
waveguide.  Assuming  that  9  -  x 12  for  this  case  (geom¬ 
etry  in  Fig.  1)  and  that  L  is  small  relative  to  the  charac- 
eristic  pulse  size  (to  preserve  detail  in  the  convolved  sig¬ 
nal),  the  fourth  beam  signal  radiated  is  directly  pro¬ 
portional  to  the  instantaneous  overlap  of  beams  1  and  2. 
Mathematically,  the  time  evolution  of  the  signal  beam  can 
be  written  as  Ut(t)  oc  {  £/,(f  -  2r)  C/:(t)  U^dr  where  it 
has  been  explicitly  assumed  that  the  control  beam  is 
"turned  on."  For  this  application,  the  input  signal  length 
cannot  be  larger  than  the  waveguide  dimension  which 
limits  the  process  to  pulses  of  100  ps  or  less.  Note  also 
the  time  compression  by  a  factor  of  two. 

Another  potential  application  of  four- wave  mixing  to 
real-time  processing  is  time  reversal  of  an  optical  wave¬ 
form,  as  illustrated  in  Fig.  3(b).  Here  9  =  0  and  beams 
1  and  2  are  the  waveform  to  be  reversed  and  a  control 
beam  of  constant  amplitude,  respectively.  Beam  3  is  a  5 
function  in  time  (very  short  pulse)  and  4  is  the  time  re¬ 
versed  signal.  A  crucial  factor  is  that  beams  1  and  3  can 
be  propagated  in  different  modes.  If  the  mode  number  of 
beam  3  is  larger  than  that  of  beam  1 ,  beam  3  can  overtake 
beam  1  from  "behind."  Thus,  beam  4  is  "read  out”  from 


Fig.  4.  Experimental  setup  used  for  the  observation  of  optical  phase  con¬ 
jugation  with  surface  plasmons.  The  inset  shows  (he  wavevectors  of  the 
three  incoming  waves  in  the  plane  of  the  interface  (BS  beam  splitters. 
VDL.  variable  delay  lines.  PA:  pnsm-silver  film  assemblv  mounted  on 
a  rotary  table.  L:  1  m  focal  length  lens  A:  0  5  mm  diameter  aperture. 
PMT:  photomultiplier  tube.  PD  photodiode  (101  l 

trailing  edge  to  leading  edge,  which  corresponds  to  time 
reversal  of  the  pulse  envelope. 

VI.  Surface  Plasmon  Experiments 
A  variety  of  experiments  have  been  earned  out  by  Nunzi 
and  Ricard  [10]  on  phase  conjugation  with  surface  plas¬ 
mons  in  attenuated  total  reflection  (ATR)  geometries.  In 
the  visible  and  near  infrared,  the  propagation  distance  for 
surface  plasmons  is  typically  tens  of  micrometers  and  the 
nonlinear  mixing  takes  place  under  a  single  coupling 
prism. 

In  the  first  set  of  experiments,  three  incident  beams  were 
coupled  into  surface  plasmon  modes  as  shown  in  Fig.  4 
and  the  conjugate  beam  was  sampled  via  a  beam  splitter. 
For  a  silver  film  on  a  glass  prism  probed  in  the  Kretsch- 
mann  geometry,  the  four-wave  mixing  signal  observed 
was  due  to  heating  of  the  metal  film.  This  was  verified  by 
pump-probe  experiments  on  the  large  spacing  nonlinearly 
produced  grating.  With  pulses  28  ps  long,  0.6  ns  decay 
times  were  measured,  in  good  agreement  with  calcula¬ 
tions.  When  liquid  ethanol  was  placed  between  the  cou¬ 
pling  prism  and  the  silver  film,  thermal  gratings  were 
formed  in  both  the  liquid  and  metal,  and  interference  ef¬ 
fects  were  observed.  The  conjugate  reflectivity  in  these 
experiments  was  on  the  order  of  10-4.  This  case  has  been 
analyzed  not  only  by  the  experimental  group  [10],  but  also 
by  Ujihara  for  nonlinearity  in  the  dielectric  medium  [22] 
and  in  metal  [23]  (modeled  as  a  free  electron  gas). 
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•  In  the  second  set  of  experiments,  two  counterpropagat- 
ing  surface  plasmons  were  excited  under  a  coupling  prism 
to  produce  a  refractive  index  grating.  This  grating  was 
then  used  to  couple  in  a  third  wave,  incident  from  the  air 
side  of  the  metal  film,  into  a  surface  plasmon  mode.  The 
nonlinearly  in-coupled  surface  plasmon  was  observed 
when  it  was  coupled  out  by  the  prism.  The  mixing  of  the 
two  counterpropagating  surface  plasmon  modes  was  also 
found  to  generate  high-frequency  acoustic  waves. 

■  There  has  been  some  speculation  that  the  efficiency  of 
DWFM  by  surface  plasmon  could  be  increased  by  using 
the  long-range  mode  guided  by  metal  films  of  thickness 
on  the  order  of  the  optical  field  skin  depth  (6],  (14).  Al¬ 
though  it  appears  theoretically  that  the  efficiencies  should 
be  larger  [6],  [14],  it  has  been  argued  by  Nunzi  and 
Ricard  (10)  that  using  such  thin  films  would  not  eliminate 
the  heating  effects  which  dominated  their  experiment. 

VII.  Summary 

The  single  factor  which  has  held  ba^  k  the  development 
of  degenerate  four-wave  mixing  in  guided  wave  structures 
has  been  the  lack  of  nonlinear  materials  suitable  for  fab¬ 
rication  in  waveguide  form.  In  the  last  few  years,  a  num¬ 
ber  of  new  material  systems  have  been  identified  which 
can  be  fabricated  into  thin  film,  and  probably  fiber  form. 
They  are  essentially  semiconductor-(24)-[27)  and  metal- 
doped  (28)  glasses  and  thin  organic  films  (29).  All  of  these 
materials  have  very  fast  relaxation  times,  with  the  semi¬ 
conductor  glasses  being  the  slowest  with  response  times 
(electronic)  of  tens  of  picoseconds.  At  the  present,  the 
semiconductor-doped  glasses  are  the  only  ones  with  suf¬ 
ficient  nonlinearity  to  result  in  waveguide  four-wave  mix¬ 
ing  with  subwatt  peak  power  levels.  However,  improve¬ 
ments  are  expected  in  the  field  of  nonlinear  organic 
materials  so  that  they  become  competitive  with  the  doped 
glass  materials. 

With  such  highly  nonlinear  and  ultrafast  materials, 
some  of  the  promising  applications  to  all-optical  signal 
processing  should  become  feasible,  and  we  look  forward 
to  these  possibilities  with  great  anticipation. 
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The  properties  of  planar  waveguides  and  directional  couplers  fabricated  by  potassium/sodium 
ion  exchange  in  semiconductor-doped  glasses  are  described. 


A  new  spectrum  of  ultrafast  (picosecond)  all-optical 
devices  becomes  possible  when  an  optical  waveguide  con¬ 
tains  a  medium  characterized  by  an  intensity-dependent  re¬ 
fractive  index. 1  Although  many  interesting  device  concepts 
have  been  proposed,'  progress  has  been  hampered  by  the 
lack  of  nonlinear  materials  with  both  large  nonlinearities 
and  picosecond  response  times,  and  from  which  optical 
waveguides  can  easily  be  fabricated.  To  date,  only  a  nonlin¬ 
ear  coherent  coupler,  first  proposed  by  Jensen,2  has  been 
demonstrated  with  cw  excitation  in  strain-induced  wave¬ 
guides  in  GaAs/AlGaAs  multiple  quantum  well  material.3'4 
This  device  consisted  of  two  parallel,  channel  waveguides 
and,  when  light  is  initially  inserted  into  one  channel,  the 
transfer  efficiency  into  the  second  channel  depends  upon  the 
initial  power  in  the  incidence  channel.  A  further  refinement 
is  a  nonlinear  Mach-Zehnder  interferometer  as  described  by 
Lattes  et  al ,5  Recently  it  has  been  shown  that  semiconduc¬ 
tor-doped  glasses  (color  filters),  which  are  glasses  coatain- 
ing  small  ( 10-100  nm)  crystallites  of  CdS,  Se,  _x  .have  use¬ 
fully  large  nonlinearities4-7  (n2  =  10“ 14  m2/W)  for  guided 
wave  devices,  coupled  with  picosecond  relaxation  times7 
(r<10_"  s ) .  For  example,  the  critical  power  Pe  foranon- 
linear  coherent  coupler2  acting  as  an  all-optical  switch 
should  be  approximately  1  W  for  a  channel  cross-sectional 
area  of  10  fim 2  and  a  linear  transfer  length  of  5  mm.  Since 
such  glasses  are  relatively  inexpensive  and  can  be  made  with 
10%  concentrations  of  sodium,  they  are  promising  candi¬ 
dates  for  waveguide  fabrication  by  ion  exchange,  one  of  the 
more  simple  waveguide  fabrication  technologies.*'9  In  this 
letter  we  describe  the  fabrication  and  the  waveguiding  char¬ 
acteristics  of  both  planar  waveguides10  and  directional  cou¬ 
plers  in  semiconductor-doped  glasses. 

The  color  filter  glasses  used  in  this  study  consist  of  a 
bulk  soda-potash-lime  glass  matrix  doped  with  small  crys¬ 
tallites  of  the  semiconductor  compound  CdSz  Se,  _  x  with  x 
chosen  to  give  a  red  colored  glass  with  a  band  edge  around 
390  nm.  The  nonlinearity  is  resonant  with  the  band  gap  and 
this  wavelength  is  convenient  for  subsequent  characteriza¬ 
tion  of  the  nonlinear  properties  of  the  devices  with  a  mode- 
locked  rhodamine  6G  dye  laser.  The  glass  content  is  9.5%  by 
weight  sodium  and  3.3%  by  weight  potassium.  The  sub¬ 
strate  index  was  measured  using  an  Abbe  refractometer  to  be 
1 .5 1 83  ±  0.0001  at  632.8  nm.  Planar  waveguides  were  fabri¬ 
cated  by  immersing  the  filter  into  a  bath  of  molten  potassium 
nitrate  at  340  *C  ( this  temperature  just  above  the  melting 
point  of  334  *C  was  found  to  give  the  most  reproducible  re¬ 


sults  )  in  which  sodium  ions  leave  the  glass  to  be  replaced  by 
potassium  ions  from  the  melt.  This  exchange  results  in  a 
region  of  higher  refractive  index  An  of  around  0.009  at  the 
surface  which  decreases  monotonically  with  depth.  To  fabri¬ 
cate  stripe  waveguides,  an  aluminium  diffusion  mask  is  de¬ 
posited  on  the  glass:  the  mask  has  apertures  through  which 
selective  ion  exchange  can  take  place  thus  defining  the  wave¬ 
guide  dimensions. 

The  fabrication  of  single-mode  planar  and  channel 
waveguides  requires  a  detailed  investigation  of  the  ion-ex¬ 
change  process.  A  characterization  study  of  K^/Na*  ion 
exchange  in  color  filter  glass  was  undertaken  similar  to  those 
reported  in  Ref.  8  (Ag+/Na+  ion  exchange)  and  Ref.  11 
(K+/Na+  ion  exchange)  for  soda-lime  microscope  slide 
glass  and  Ref.  12  (K^/Na+  ion  exchange)  for  BK-7  glass . 
A  number  of  planar  waveguide  samples  were  prepared  in 
which  the  diffusion  times  t  were  varied  between  the  range  1 3 
min  to  24  h.  For  every  sample  the  effective  index  of  each 
mode  was  determined  by  the  prism  coupling  method  to  an 
accuracy  of  ±  0.0001  with  a  He-Ne  laser  operating  at  wave¬ 
length  of  632.8  nm.  The  results  for  both  TE  and  TM  polari¬ 
zations  are  shown  in  Fig.  1. 

The  theoretical  dispersion  curves  which  produced  an 
excellent  fit  to  the  data  were  derived  as  follows.  First,  the 
inverse  Wentzel-Kramers-Brillouin  (WKB)  method13  was 
used  to  obtain  the  mode  turning  points  and  also  the  surface 
refractive  index  n,,  consistently  found  to  be  1.5272  for  TE 
modes  and  1.3305  for  TM  modes.  In  order  to  find  depth 
index  profiles  which  accurately  fit  these  points,  the  finite 
difference  method  was  used  to  calculate  solutions  to  the  non¬ 
linear  diffusion  equation  in  one  dimension  given  by 


dc  _  d  f 

D ,  dc\ 

(la) 

dt  dx  \ 

1  —ac  dx)' 

C  =  Ci/Cq, 

a  -  1  -  (Z>,/2>0)  , 

(lb) 

where  c  is  the  normalized  potassium  ion  concentration  (c, 
being  the  concentration  of  the  potassium  ions  and  c0  the 
total  number  of  diffusing  potassium  and  sodium  ions  in  the 
glass)  and  D ,  and  D0  are  the  self-diffusion  constants  of  the 
potassium  and  sodium  ions,  respectively.  The  surface  con¬ 
centration  ratio  was  determined  by  laser  induced  ion  mass 
analysis  (LIMA)  of  the  exchanged  glass  surface  layer  to  be 
9:1  K  +  /Na*,  giving  a  boundary  condition  of  c  =  0.9.  The 
refractive  index  is  assumed  to  vary  linearly  with  the  potas¬ 
sium  ion  concentration.  The  shape  and  depth  of  the  index 
profiles  are  adjusted  by  varying  the  parameters  a  and  D , : 
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FIG-  I.  Theoretical  dispersion  curves  calculated  using  the  Wentzet- 
Kramers-Bnllouin  ( WKfl)  method  for  the  optimum  index  profiles  com* 
pared  to  measured  mode  effective  index  values  for  samples  fabricated  at 
T  =»  340  ‘C.  (a)  TE  modes  and  (b)  TM  modes. 


values  of  a  =  0.3  and  =  4.95x10*“  m2s~‘  were 
found  to  give  the  best  least-squares  fit  for  the  TE  modes, 
while  a  =  0.3  and  D ,  =  4.50X  10  *  “  m2  s  '  1  give  the  best 
fit  for  the  TM  modes.  The  apparent  anisotropy  of  ion-ex¬ 
changed  waveguides,  as  suggested  by  the  discrepancy 
between  the  diffusion  constants  and  the  surface  index  values 
for  the  TE  ana  TM  modes,  is  thought  to  be  due  to  stresses 
introduced  into  the  glass  by  the  different  sizes  of  the  two 
ionic  species  being  exchanged.13 

Using  the  best  fit  refractive  index  profiles  and  defining 
d ,  to  be  the  depth  where  the  index  change  is  equal  to  \/e  of 
the  surface  index  change,  a  linear  relationship  between  d, 
and  t 1/2  was  obtained.  From  the  slope,  the  “effective”  diffu¬ 
sion  constant  [defined  by  the  relation  (d,  =  (Dft)w2  ]  is 
given  by  D,  =  10.80X  10" 16  m2  s~‘  for  TE  modes  and 
D,  =9.32x10*  16  m2  s*1  for  TM  modes.  Finally,  the 
WKB  method  was  employed  using  the  diffusion  parameters 
obtained  above  to  calculate  the  curves  fitted  to  the  data  in 
Fig.  1 .  The  excellent  fit  indicates  that  the  refractive  index 
profiles  and  the  Dt  values  determined  are  reasonably  accu¬ 
rate.  This  agreement  is  particularly  important  for  predicting 
the  short  diffusion  times  required  for  single-mode  wave¬ 
guides.  A  diffusion  time  of  1  h  was  chosen  from  these  results 
as  the  optimum  time  for  the  fabrication  of  single-mode  wave¬ 
guides. 
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FIG.  2.  Schematic  diagram  of  directional  coupler  pattern  used  to  fabricate 
devices.  This  pattern  was  repeated  three  times  on  the  mask  for  different 
waveguide  separations  (d)  of  (a)  4.0,  (b)  4.3.  and  (c)  5.0  /am.  The  interac¬ 
tion  length  was  varied  in  steps  of  0.3  mm  over  the  ranges  of  (a)  1.0-3.00 
mm,  (b)  2.0-4.0  mm,  and  (c)  3.0-5.Q  mm.  The  waveguide  width  (w)  was 
nominally  4.0  ftm  for  ail  the  couplers  on  the  mask. 


The  theoretical  effective  index  values  of  stripe  wave¬ 
guides  and  the  transfer  length  for  directional  couplers  were 
calculated  using  the  numerical  model  derived  by  Walker 
et  al .’  First,  the  model  is  used  to  solve  the  nonlinear  diffu¬ 
sion  equation  in  two  dimensions  to  give  the  refractive  index 
profile  rt(xy),  then  the  Helmholtz  wave  equation  is  solved 
for  such  a  region  using  the  variational  metnod  ( the  model  is 
restricted  to  quasi-TE  modes)  to  give  the  mode  effective  in¬ 
dex  values. 

For  a  diffusion  time  of  1  h  at  340  *C  the  slab  effective 
index  is  1.5205.  The  theory  predicts  that  stripe  waveguides 
formed  with  the  same  diffusion  parameters  will  be  single 
mode  for  aperture  widths  w  of  3-5  /im  with  effective  index 
values  from  1.5189  to  1.5196.  The  mask  design  used  to  ex¬ 
perimentally  verify  directional  coupler  calculations  is  shown 
in  Fig.  2. 

Several  samples  were  fabricated  using  this  mask  and  the 
edges  polished  to  facilitate  end-fire  coupling  ( this  process 
removed  0. 5  mm  from  the  interaction  length  of  each  device ) . 
Light  from  a  He-Ne  laser  was  coupled  in  turn  into  each  de- 
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FIG.  3.  Plot  of  the  measured  normalized  coupled  power  <  I  —  P  P , )  vs 
interaction  length  (£, )  for  the  three  sets  of  similar  couplers  on  a  sample. 
The  best  fit  cos*  ( trr/2ir )  curves  are  also  shown.  The  results  are  shown  in 
Table  I.  (The  lengths  are  0.5  mm  less  than  the  lengths  in  Fig.  2  because  this 
amount  is  lost  in  the  polishing  process. ) 

Cutlon  mat.  1404 


J 


TABLE  l.  Comparison  of  theoretical  and  measured  coupler  transfer  length 
for  devices  prepared  at  r  =  340  ‘C  for  1  h. 


Coupler 

set 

Waveguide 
width  w 

(pm) 

Waveguide 
separation  d 
(Aim) 

Transfer  length  L,  ( mm ) 
measured  theoretical 

a 

4.2 

4.2 

3.5  ±0.25 

3.5 

b 

4.2 

4.6 

5.2  ±  0.25 

5.5 

c 

4.2 

5.3 

7.7-0.35 

7.9 

All  guide  widths  and  separations  measured  to  an  accuracy  of  ±0.1  /im. 


vice  on  a  sample  at  the  single-guide  end  and  the  power  in 
each  of  the  output  guides  was  determined  by  focusing  the 
near-field  image  onto  the  camera  of  a  calibrated  vidicon  sys¬ 
tem.  Figure  3  shows  the  power  remaining  in  the  input  guide 
versus  the  device  interaction  length  for  the  three  sets  of  de¬ 
vices.  The  best  fit  cos2  ( rrz/2Lc )  curves  are  also  shown  to 
estimate  the  coupling  length.  Table  I  gives  the  actual  device 
parameters  and  measured  transfer  lengths  compared  with 
the  predicted  transfer  lengths  calculated  from  the  theoretical 
model.  Very  good  agreement  was  obtained. 

In  summary,  good  quality  single-mode  directional  cou¬ 
plers  can  be  made  by  K  *  /Na  *"  ion  exchange  in  semiconduc¬ 
tor-doped  glasses.  The  waveguiding  properties  of  a  direc¬ 
tional  coupler  can  also  be  predicted  with  the  model 
presented  here.  Experiments  are  currently  under  way  to 
measure  the  nonlinear  properties  of  the  directional  couplers, 
as  well  as  the  Mach-Zehnder  devices. 

This  research  was  supported  in  the  United  Kingdom  by 
the  Science  and  Engineering  Research  Council  and  in  the 


United  States  by  the  National  Science  Foundation  (ECS- 
88501249),  the  Army  Research  Office  (DAAG29-85-K- 
0173 ),  NSA  via  the  Joint  Services  Optics  Program,  and  the 
Air  Force  Office  of  Scientific  Research  ( AFOSR-84-0277). 
The  laser  induced  ion  mass  analysis  ( LIMA )  measurements 
were  performed  courtesy  of  Plessey  Ltd. 


‘Recently  reviewed  in  G.  1.  Stegeman  and  C  T  Seaton.  Appl  Phys.  Rev 
( in  J.  Appl.  Phys. )  58.  R57  ( 1985). 

:S.  M.  Jensen,  IEEE  J  Quantum  Electron.  QE-18.  1580  ( 1982). 

>p.  Li  ICam  Wa.  J.  E.  Such.  N  J  Mason.  J  S  Roberts,  and  P  N.  Robson. 
Electron.  Lett.  21.26(1985) 

4P.  Li  JCam  Wa.  J  H  Marsh.  P  N  Robson.  J  S  Roberts,  and  N.  J.  Mason. 
Proceedings  of  the  2nd  Conference  on  Integrated  Optical  Engineering. 
Cambridge.  SP1E  (The  Society  of  Photo-Optical  Instrumentation  Engi¬ 
neers),  edited  by  S.  Siram,  578.  1 10  ( 1985). 

5A.  Lattes.  H.  A.  Ha  us.  F  J  Leonberger.  and  E  P  Ippen.  IEEE  J  Quan¬ 
tum  Electron.  QE-19.  1718  (1983) 

*R.  K.  Jam  and  R.  C.  Lind,  J  Opt.  Soc.  Am.  73.  647  ( 1983) 

7S.  S.  Yao.  C.  ICaraguleff.  A.  Gabel.  R.  Fortenberry.  C.  T.  Seaton,  and  G  1 
Stegeman.  Appl  Phvs.  Lett.  46,  801  ( 1985). 

"G.  Stewart,  C  A.  Millar.  P  J  R  Layboum,  C.  D  W  Wilkinson,  and  R 
M.  De  La  Rue,  IEEE  J.  Quantum  Electron  QE-13.  192  ( 1978). 

H.  G.  Walker.  C.  D  W  Wilkinson,  and  J  .  A  .  H  Wilkinson,  Appl  Opt  22 
1923  (1983). 

loC.  N.  Ironside.  J  F  Duffy.  R  H  Hutchins.  W  C.  Banyai.  C  T  Seaton, 
and  G.  I.  Stegeman,  Technical  Digest  of  the  5th  International  Conference 
on  Integrated  Optics  and  Optical  Communications  and  of  the  1  Ith  Euro¬ 
pean  Conference  on  Optical  Communications,  Venice.  Italy.  1985,  p.  237 
"G.  L.  Yip  and  J.  Albert.  Opt.  Lett  10,  151  (1985) 

1JJ.  E.  Gortych  and  D.  G.  Hall.  Opt.  Lett.  11  100  ( 1986) 
l5J.  Finak.  H.  Jerominek.  Z  Opilski.  and  K  Wouala.  Opt  Appl  XII.  1 1 
(1982). 


1405 


Appl  Phy*.  Lott..  Vol  49.  No  21  24  November  1986 


Cuiien  et  ai 


1405 


0  54  8 


AFOSR.TK-  8  7- 

Gaussian  beam  excitation  of  TE0  nonlinear  guided  waves 

E.  M.  Wright,  G.  I.  Stegeman,  and  C.  T.  Seaton 

Optical  Sciences  Center  and  Arizona  Research  Laboratories.  University  of  Arizona.  Tucson.  Arizona  85721 

J.  V.  Moloney 

Department  of  Physics.  Heriot-Watt  University.  Ricarton.  Edinburgh.  EH  14  4AS.  United  Kingdom 
(Received  16  May  1986;  accepted  for  publication  1  July  1986) 

Selective,  efficient  excitation  of  TE„  nonlinear  guided  waves  is  demonstrated  numerically.  For 
a  thin  film  bounded  by  two  self-focusing  media,  three  different  field  distributions 
corresponding  to  the  same  flux  level  can  be  excited  independently  by  suitable  Gaussian  input 
beams. 


To  date  several  integrated  optics  devices  have  been  pro¬ 
posed  utilizing  nonlinear  waveguides  ( NLWG’s ) ,  for  exam¬ 
ple,  upper  and  lower  threshold  devices,  nonlinear  couplers, 
optical  switches,  etc. 1  Such  proposals  have  invariably  been 
based  on  the  NLGW  dispersion  curves  (guided  wave  flux 
versus  effective  index)2  even  in  the  absence  of  data  on  the 
stability  of  the  various  branches  and  of  the  possibility  of  ex¬ 
ternally  exciting  the  NLGW’s.  Moloney  et  al.3A  and  Jones 
and  Moloney3  have  recently  addressed  the  stability  question 
both  numerically  and  analytically,  and  a  comprehensive  pic¬ 
ture  is  now  emerging.  In  contrast,  the  problem  of  external 
excitation  of  NLGW’s  (in  the  waveguide  context)  has  not 
been  discussed  so  far,  and  is  the  subject  of  this  letter.  Specifi¬ 
cally,  we  show  that,  for  a  fixed  input  flux,  selective  excitation 
of  the  corresponding  TE^  NLGW’s  by  end  firing  of  appro¬ 
priately  tailored  Gaussian  beams  is  indeed  possible.  Implica¬ 
tions  for  device  design  concepts  are  also  discussed. 

We  treat  the  case  of  TE  waves  of  frequency  to  in  a  planar 
NLGW  as  shown  in  Fig.  1.  For  this  preliminary  study  the 
electric  field  is  assumed  homogeneous  in  the  y  direction. 
Then  writing  the  electric  field  as 

E(r,t)=\y[W(xj)e“B'-M' +c.c.)  (1) 

gives  in  the  usual  slowly  varying  envelope  approximation, 

[02 -n2(x,\W\2))W=Q,  (2) 

02  ax' 

where  the  refractive  index  in  the  various  media  is  described 
by 

n2{x.\iV\2)  =n2r  +  ar\W(xj)\2.  r  =  c,/4,  (3) 

and  the  subscripts  c,  f ,  s  refer  to  the  cladding  {x<.  —  d), 
film  (  —  d<.x<d),  and  substrate  (x>d),  respectively.  Here 
the  coordinates  x  and  z,  and  d  are  in  units  of  k  1  =  c/a.  The 
NLGW’s  of  the  system  are  solutions  of  Eq.  (2)  such  that 
W(xs)  =  fV0(x).  Equation  (2)  along  with  the  boundary 
conditions  of  continuity  of  W  and  dW  /dx  at  the  interfaces 
(jc  =  ±d)  can  then  be  used  to  generate  the  dispersion 
curves  for  the  system.  Note  that  due  to  the  nonlinearity  as¬ 
sumed  in  this  problem,  several  NLGW  solutions  may  exist 
for  a  given  1,2  An  example  of  the  dispersion  curves  for  the 
TE,,  NLGW’s  of  a  waveguide  with  nonlinear  cladding  and 
substrate,  and  a  linear  film,  is  shown  in  Fig.  2  ( we  use  the 
parameter  values  given  in  this  figure  throughout).  We 
choose  self-focusing  nonlineanties  since  this  case  contains 
by  far  the  richest  set  of  new  phenomena. 12  The  regions 
marked  by  the  dashed  lines  are  unstable  under  propagation. 


and  the  remaining  regions  are  stable  and  are  labeled  I,  II,  and 
III,  respectively.'  The  inset  beside  each  stable  region  indi¬ 
cates  the  field  distribution  of  the  NLGW  on  that  branch.  For 
example,  in  region  II  the  NLGW  is  associated  with  self-fo¬ 
cusing  in  the  medium  with  the  lower  nonlinear  coefficient, 
whereas  in  III  it  is  associated  mainly  with  self-focusing  in  the 
other  bounding  medium.  At  high  flux  levels  these  solutions 
degenerate  into  nonlinear  surface  polaritons  at  their  respec¬ 
tive  interfaces,  each  independent  of  the  other  interface. 

To  investigate  the  excitation  of  NLGW’s  we  have  solved 
Eq.  (2)  with  initial  data  corresponding  to  a  collimated 
Gaussian  input  beam; 

Wm(x)  =  WNe-'x-‘°)‘/'*,  (4) 

where  WN  is  used  to  adjust  the  input  beam  flux,  xQ  is  the 
displacement  of  the  beam  center  from  the  guide  center,  and 
w0  is  the  Gaussian  spot  size.  This  boundary  condition  corre¬ 
sponds  to  end  firing  a  Gaussian  beam  onto  the  NLGW.  Note 
that  we  ignore  any  nonlinear  reflection  that  may  occur  at  the 
air-waveguide  interface.  This  is  valid  since  the  input  beam  is 
incident  normally  onto  the  interface  whereas  nonlinear  ef¬ 
fects  occur  only  for  angles  close  to  the  critical  angle  for  total 
internal  reflection.6  7 

For  our  numerical  experiment  we  have  fixed  the  input 
flux  (S,n  =  0. 165),  and  as  indicated  by  the  horizontal  line  in 
Fig  2,  there  are  three  distinct,  stable  NLGW  solutions.  For 
each  solution  we  examined  the  NLGW  profile  and  estimated 
values  for  x0  and  w0  such  that  the  input  Gaussian  would 
resemble  the  corresponding  NLGW.  For  example,  for  the 


FIG  1  Nonlinear  waveguide  geometry,  a  linear  thin  film  of  thickness  Id  is 
bounded  by  two  self-focusing  nonlinear  media 
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FIG.  2.  Guided  wave  flux  (S)  vs  effective  index  (  P)  for  two  self-focusing 
bounding  media.  The  solid  and  dashed  lines  indicate  stable  and  unstable 
regions,  respectively,  and  the  insets  show  the  nature  of  the  NLGW’s  on  the 
various  branches.  Parameter  values  are  nc  —  n,  »  1.53,  nr  m  1.57, 
a,  »  10-’,  a,  -  8x  10-’,  a,  -  0,  andd  -  8. 

solution  on  branch  I  of  Fig.  2  we  chose  iv0  =  10.0  and 
x„  =  0.0,  whereas  for  branch  II  we  have  w0  =  6.0  and 
x0  =  10.3.  The  evolution  of  the  Gaussians  for  branches  I— III 
are  shown  in  Figs.  3(a)-3(c),  respectively.  In  all  cases  the 
Gaussians  evolve  toward  the  NLGW  they  were  tailored  to 
excite.  This  was  verified  by  detailed  examination  of  the  final 
beam  profile  and  the  corresponding  NLGW.  We  have  there¬ 
fore  demonstrated  that  selective  excitation  of  the  TEq 
NLGW’s  at  a  fixed  flux  level  using  a  Gaussian  beam  is  possi¬ 
ble.  These  results  indirectly  imply  that  selective  excitation  is 


id)  x 


FIG.  3.  Evolution  of  Gaussian  input  beams  for Sm  >  0. 165  and  (a)  x„  =  0, 
ui0  »  10.0.  (b)  x0  «  10.3,  ui„  —  6.0;  (c)  xn  x  —10.6,  u;0  »  4.0,  (d) 
x„  »  -  9.0,  u/„  a  4.0.  The  propagation  coordinate  is  in  units  of  free-space 
wavelengths,  and  the  vertical  lines  are  the  waveguide  boundaries. 


also  possible  if  only  the  cladding  is  nonlinear,  in  that  case  the 
only  difference  is  that  branch  II  vanishes  and  there  are  only 
two  NLGW’s.  We  have  furthermore  verified  numerically 
that  Gaussian  beam  excitation  is  possible  for  many  values  of 
the  ratio  ac/a,.  The  results  presented  here  are  representa¬ 
tive  of  our  general  findings. 

A  new  question  immediately  arises  from  these  results, 
namely,  how  sensitive  is  the  coupling  to  the  initial  Gaussian 
beam  profile?  First,  we  stress  that  our  estimates  for  z0  and  u>0 
came  from  visual  examination  of  the  NLGW's  which  indi¬ 
cates  that  there  is  no  exceptional  sensitivity  to  initial  condi¬ 
tions.  However,  if  the  input  beam  profile  is  too  different  from 
a  NLGW  it  decays  mainly  via  emission  of  radiation  and/or  a 
soliton  which  is  ejected  into  the  nonlinear  medium  with  the 
larger  nonlinear  coefficient.4  (If  both  media  are  identical, 
two  solitons  are  emitted,  one  into  each  medium. )  An  exam¬ 
ple  of  the  soliton  decay  is  shown  in  Fig.  3(d)  which  corre¬ 
sponds  to  the  situation  in  Fig.  3(c)  but  with  the  beam  center 
misaligned  from  x0  =  —  10.6  tox0  =  —  9.0.  More  details  of 
these  decay  processes  will  be  published  elsewhere. 

These  results  show  that  the  input  beam  profile  is  a  signif¬ 
icant  control  parameter  in  considerations  of  NLGW  excita¬ 
tion.  The  input  flux  alone  is  not  sufficient  to  characterize  the 
input  This  has  implications  for  device  design  concepts.  For 
example,  previous  proposals  for  optical  power  limiting  using 
NLGW's  have  relied  on  the  fact  that  for  certain  configura¬ 
tions  (e.g.,  a  waveguide  with  a  self-defocusing  cladding 
which  starts  above  cutoff),  there  exist  no  NLGW’s  above  a 
certain  critical  flux.1  For  the  case  shown  in  Fig.  2,  we  can 
couple  a  Gaussian  beam  (x0  =  0.0,  w0  =  10.0)  into  branch  I 
very  efficiently.  (The  system  is  essentially  linear  on  this 
branch. )  However,  for  input  fluxes  above  which  this  branch 
does  not  exist,  our  input  beam  does  not  match  the  solutions 
on  the  two  remaining  branches  (II  and  III)  and  the  input 
decays  via  emission  of  radiation  and/or  solitons,  leaving  a 
NLGW  on  branch  I  of  reduced  flux.  This  process  produces 
threshold  behavior  and  optical  limiting.  Without  proper 
consideration  of  the  input  profile,  this  action  cannot  be  un¬ 
derstood  since  NLGW’s  exist  for  all  possible  input  fluxes. 

In  conclusion,  we  have  demonstrated  that  selective  exci¬ 
tation  of  stable  TEq  NLGW’s  using  a  Gaussian  input  beam  is 
possible.  Indeed  it  seems  that  simply  matching  the  input 
beam  to  a  particular  NLGW  is  sufficient  to  guarantee  selec¬ 
tive  excitation. 
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We  demonstrate  numerically  that  external  excitation  of  a  nonlinear  waveguide  can  produce 
sequential  threshold  behavior  via  multisoliton  emission  from  the  waveguide.  This  behavior  is  simi¬ 
lar  to  that  predicted  to  occur  at  a  nonlinear  interface. 


Following  the  pioneering  work  of  Kaplan1  on  the  re¬ 
flection  of  a  plane  wave  from  a  linear-nonlinear  interface, 
several  papers  appeared  on  this  subject. 2  This  interest 
stemmed  from  the  possibility  of  bistable  reflection  from 
such  an  interface. 1  Although  it  is  now  generally  accepted 
that  bistable  reflection  is  not  possible  (as  opposed  to  a 
hysteretic  response)  several  of  the  results  from  the  theoret¬ 
ical  studies,  such  as  the  nonlinear  Goos-Hiinchen7  g  effect 
and  the  transmission  of  self-focused  channels  (or  solitons) 
through  the  interface,2  7  are  of  interest  in  their  own  right. 
To  the  best  of  our  knowledge  only  the  former  effect  has 
received  further  attention. 

In  this  Brief  Report  we  consider  the  external  excitation 
of  a  nonlinear  waveguide  (NLWG)  where  the  film  and 
substrate  are  linear  but  the  cladding  displays  a  nonlinear 
refractive  index  (optical  Kerr  effect).  Intuitively  one  may 
expect  that  effects  reminiscent  of  the  nonlinear  interface 
problem  can  arise  since  the  NLWG  comprises  at  least  one 
such  interface.  Here  we  report  results  showing  the 
transmission  of  solitons  through  the  film-cladding  inter¬ 
face,  the  number  and  angle  of  emission  of  which  depends 
on  the  input  flux.  [The  input  beam  profile  is  held  fixed  as 
the  zeroth-order  transverse-electric  (TE)  mode  of  the 
linear  guide  ]  The  flux  trapped  in  the  waveguide  then 
shows  sequential  threshold  behavior  as  a  function  of  the 
input  flux  as  a  result  of  this  multisoliton  emission.  An  in¬ 
tuitive  explanation  for  the  thresholds  utilizing  the  NLWG 
dispersion  curve  for  the  system  is  given.  We  remark  that 
the  emission  of  a  soliton  from  a  NLWG  has  previously 
been  reported  as  a  route  by  which  unstable  nonlinear  guid¬ 
ed  waves  (NLGW’s)  decay. 10  However,  in  that  case  the 
input  flux  fixes  the  input  beam  profile  whereas  here  we 
hold  the  input  profile  fixed  which  is  more  in  line  with  ex¬ 
perimental  procedures.  Potential  applications  of  this  ef¬ 
fect  include  optical  limiting,  coupling  of  adjacent 
waveguides  in  the  absence  of  evanescent  field  overlap,  and 
a  light-driven  angular  scanning  element. 

We  consider  TE  waves  of  frequency  >,>  in  a  slab 
waveguide,  then  assuming  that  the  electric  field  is  homo¬ 
geneous  in  the  y  direction  < z  and  x  being  the  propagation 
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and  transverse  coordinates  in  units  of  c /to)  and  writing 
the  field  as 

E(r,r)=  yy[  tVix.zie'^  ~‘,"'  +  c.c.  ]  (1) 

yields  the  usual  slowly  varying  envelope  equation  for 
W[  x,z):u 

2i(3~  +  — ~  -[02-n:U.  i  W  :)]R'=0  .  (2) 

dz  0X" 

Here  the  refractive  index  in  the  various  media  is  given  by 

/t2(x,  W  2)  =  nY  +aY  IFi.x.z)  2,  y=c,f,s  (3) 

the  subscripts  c.f.s  referring  to  the  cladding  < x  <  -</), 
film  (—d<x<d),  and  substrate  ix>d),  respectively. 
For  the  results  presented  here  we  specifically  set 
nf=rtj  =  1.55,  tty  =1.57,  nr  =  10  :,  and  «y=a,--0, 
which  corresponds  to  a  self-focusing  Kerr-type  nonlinear 
cladding. 

The  NLGW's  are  found  as  the  solutions  of  Eq.  (2)  such 
that  IV(x.z)-  H',„(.<):12  Equation  (2)  along  with  the  con¬ 
tinuity  conditions  on  the  tangential  components  of  the 
electric  and  magnetic  fields  can  then  be  used  to  generate 
the  dispersion  curves  (guided  wave  flux  S  versus  effective 
index  13)  for  the  system.  In  Fig.  1(a)  we  show  the  disper¬ 
sion  curve  for  the  TE<,  NLGW's  of  a  NLWG  described  by 
the  parameters  given  above  (note  that  since  there  is  only 
one  nonlinear  medium  the  NLGW  is  uniquely  specified 
by  (3).  For  this  particular  configuration  the  negative-slope 
region  of  the  dispersion  curve  is  unstable  under  propaga¬ 
tion  and  is  indicated  by  a  dashed  line,  the  regions  marked 
I  and  II  are  stable. 11  The  inset  beside  each  branch  shows 
the  general  nature  of  the  solution  on  that  branch  (e.g..  on 
branch  II  the  solution  is  concentrated  around  the  film- 
cladding  interface1'),  and  the  critical  flux  Y(  is  that  lluv 
above  which  the  NLGW's  on  branch  l  cease  to  exist  for 
the  example  here  .V,  ~  0  ll>5  >. 

Our  numerical  experiment  consisted  of  launching  an  m 
pul  beam  of  fixed  profile  hut  variable  flux  omo  the 
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Some  features  of  this  soliton  emission  draw  attention. 
Firstly  and  perhaps  most  importantly,  the  first  threshold 
provides  a  highly  desirable  optical  limiter  characteristic: 
the  input-output  transfer  is  essentially  100%  right  up  to 
this  point  at  which  the  output  drops  drastically.  Presum¬ 
ably  the  switch  contrast  shown  in  Fig.  Kb)  could  be  fur¬ 
ther  increased  by  varying  the  guide  parameters.  Further¬ 
more,  the  output,  which  is  always  a  branch-I  NLGW  and 
closely  resembles  the  linear  TE„  mode  of  the  waveguide,  is 
collimated.  Thus,  in  contrast  to  previous  proposals  for 
optical  limiters  based  on  nonlinear  refractive  effects,15-17 
the  present  system  does  not  impose  any  undesirable 
phase-curvature  on  the  output  beam.  Secondly,  the  angle 
of  emission  of  the  solitons  is  a  function  of  the  input  flux. 
This  effect  was  also  observed  by  Tomlinson  el  at.  in 
their  study  of  nonlinear  interfaces.  As  pointed  out  by 
these  authors,  this  effect  has  potential  use  as  a  light- 
driven  angular  scanning  element.  Finally,  we  have  found 
that  a  soliton  emitted  from  one  NLWG  can  be  at  least 


'A.  E.  Kaplan,  Pis'ma  Zh.  Eksp.  Teor.  Fiz..  24,  132  (1976) 
(JETP  Lett.  24.  114  (1976)];  Zh.  Eksp.  Teor.  Fiz.  72,  1710 
( 1977)  [Sov.  Phys— JETP  45,  896  (1977)]. 

-N.  N.  Rosanov,  Opt.  Spektrosk.  47,  606  (1979)  [Opt.  Spectrosc. 
(LSSRi  47,  335  (1979)]. 

'P.  W.  Smith,  J.  P.  Hermann.  W.  J.  Tomlinson,  and  P.  J.  Malo¬ 
ney,  Appl.  Phys.  Lett.  35,  846  (1979). 

4D.  Marcuse,  Appl.  Opt.  19.  3130  (1980). 

'W.  J  Tomlinson,  Opt  Lett.  5.  323  (1980). 
hP.  W.  Smith.  W  J  Tomlinson,  P.  J  Maloney,  and  J.  P.  Her¬ 
mann.  IEEE  J  Ouantum  Electron.  QE-17.  340  (1981). 

W.  J.  Tomlinson,  J.  P  Gordon.  P.  W.  Smith,  and  A.  E.  Kap¬ 
lan,  App.  Opt.  21.  2041  (1982). 

'P  \V  Smith  and  \V.  J  Tomlinson,  IEEE  J.  Ouantum  Electron. 
OE-20.  30  (1984  . 

“N.  N.  Akhemdiev.  V.  I  Korneev,  and  Yu.  V.  Kuz'menko,  Zh. 
Eksp  Teor.  Fiz.  88.  107  (1985)  [Sov.  Phys— JETP  61.  62 
■  |9X5i. 

lnJ  V.  Moloney,  J  Ariyasu.  C.  T.  Seaton,  and  G.  I.  Stegeman, 


partially  trapped  by  a  second  adjacent  guide.  This  offers 
the  possiblity  of  directional-coupler  action  without  the 
need  for  evanescent  field  overlap  of  the  guides.  Further¬ 
more,  since  soliton  emission  occurs  only  above  a  given 
threshold  flux,  coupling  between  the  guides  will  experi¬ 
ence  the  same  threshold  behavior  which  could  be  used  to 
construct  logic  elements. 

In  conclusion,  we  have  shown  that  external  excitation 
of  a  NLWG  can  produce  sequential  threshold  behavior 
through  multisoliton  emission,  and  a  simple  explanation 
for  the  threshold  behavior  has  been  advanced  based  on  the 
NLGW  dispersion  curves. 
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